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ABSTRACT 

This paper provides a complete theoretical treatment of the point-mass lens perturbed by 
constant external shear, often called the Chang-Refsdal lens. We show that simple invariants 
exist for the products of the (complex) positions of the four images, as well as moment sums of 
their signed magnifications. The image topographies and equations of the caustics and critical 
curves are also studied. We derive the fully analytic expressions for pre-caustics, which are 
the loci of non-critical points that map to the caustics under the lens mapping. They constitute 
boundaries of the region in the image domain that maps onto the interior of the caustics. The 
areas under the critical curves, caustics and pre-caustics are all evaluated, which enables us 
to calculate the mean magnification of the source within the caustics. Additionally, the exact 
analytic expression for the magnification distribution for the source in the triangular caustics is 
derived, as well as a useful approximate expression. Finally, we find that the Chang-Refsdal 
lens with additional convergence greater than unity (the 'over-focusing case') can exhibit 
third-order critical behaviour, if the 'reduced shear' is exactly equal to V3/2, and that the 
number of images for N point masses with non-zero constant shear cannot be greater than 
5N-1. 



1 INTRODUCTION 

The so-called Chang-Refsdal lens was put forward to describe the 
lensing effe cts of stars in a background galaxy <Chang & Refsdall 
1979, 1984). The star acts as a point-mass lens, while the galaxy 
provides a background perturbation field, which can be approx- 
imated by a tidal term to lowest order. The Chang-Refsdal lens 
is therefore the simplest description of lensing by a tidally per- 
turbed point-mass. As such, it has found widespread astronom- 
ical applicat ions - e.g. in the m odelling of binary micr olensing 
lightcurves iGould & Loeblll992t iGaudi & Gould 1 19971) . in the 
analysis of the statistics of high-magnification events caused by 
stars in foreground galaxies (Schneider 1983), an d m the study of 
microlensi ng of stars arou nd the black hole in the Galactic Centre 
(Alexander & Loeb 2001), and so on. 

The Chang-Refsdal lens is of outstanding physical impor- 
tance, but it also has a number of elegant mathematical proper- 
ties, primarily because the 'deflection function' is a rational func- 
tion of position in the lens plane. Although there has already been 
much work done on the structure of the caustics a nd critical curves 
JChang & Refsdall 984 ISubramanian. Chitre. & NarasimhJl985l 
Schneider. Ehler s!& Falcdll992h . many of the properties of the 
Chang-Refsdal lens do not appear to be in the literature. The pur- 
pose of this paper is to give a compendium of those results. 

The paper is organized as follows. Section [2] presents the 
Chang-Refsdal lens utilizing the complex notation popularized by 
IWittH 199ft see also lBourassa & Kantowskill975li . In section|3| the 
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lens equation is converted into the imaging polynomial, the roots 
of which include the image locations. Provided that the source lies 
within the caustic - so that the number of images is four - then there 
exist invariants, such as the products of the image positions or the 
sums of the magnifications weighted by powers of the positions. 
Section |4] considers two cases when the lens equation is solvable 
via elementary means, namely when the source lies on one of the 
axes of the symmetry. In section|5| the pre-caustics of the Chang- 
Refsdal lens are isolated, together with the critical curves and the 
caustics. This enables us to calculate in section [6] the mean mag- 
nification of the 4-image configurations, and the exact form for a 
certain conditional distribution of magnifications. Finally, sections 
and [8] consider two generalizations of the Chang-Refsdal lens, 
namely the cases of a convergent background and of N point masses 
with shear. We find that, with the convergence greater than unity, 
the number of images is either nil or two if the 'reduced shear' is 
smaller than V3/2, nil, two, or four if it is in between V3/2 and 
the unity, and two or four if it is greater than unity. As for N point 
masses with shear, we establish that the maximum number of pos- 
sible images are bounded by 5iV - 1. 



2 THE LENS EQUATION 

The lowest order effect of the tidal field caused by the external 
mass distribution on gravitational lensing is usually described by 
a quadratic function that approximates the potential caused by the 
external masses. Here, we assume that there is no external mass 
locally, i.e., V 2 i/r ext = (zero local convergence; see section^for 
the discussion regarding the generalization when V 2 i/r cxt = 2k > 0). 
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Then, the lensing potential i^ cxt due to the external mass distribution 
may be approximated to be (frext(*i> %i) = 7r(^i - *f)/2 + y\X\X2, 
which is actually the most general form of a quadratic function 
satisfying Laplace's equation. (The constant term in the potential 
has no physical consequence whilst the linear terms lead to a con- 
stant deflection, which can be ignored by introducing the 'offset' 
between the coordinate origins for the source and the image po- 
sition.) If a point-mass lens is present under the influence of this 
tidal field, the whole system is described by the total potential 
if/ = (1/2) ln(x 2 + x^) + if/ ext , which is a superposition of the point- 
mass potential and the external potential. Here, the coordinate ori- 
gin in the lens/image plane is defined by the line of sight towards 
the lens, and the unit of the angular measurement is given by the 
'Einstein ring' radius corresponding to the mass of the point-mass 
lens. 

Then, since the deflection angle is the gradient of the lensing 
potential, we find the lens equation (y = x — Vi/0: 
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(1) 



or more compactly, utilizing the com plex number notation 
iBourassa & Kantowsldl975UWitl 19901) 



shear' (i.e., the —T 1 term) due to the point mass and the constant 
external shear y. 

We note that it is always possible to choose the direction of 
the real axis such that y is a positive real number. In particular, if 
y = lyle 2 "* 1, in a given coordinate system, the rotation of the coordi- 
nate axis by an angle <f> y with respect to the origin achieves this. The 
direction defined by is in fact one of the eigendirections of the 
shear tensor, which we shall refer to as the direction of the shear. 
The eigenvalue of the shear tensor associated with the direction of 
the shear is \y\, which is sometimes just referred to as the shear by 
itself. Since the shear tensor is real symmetric and traceless, the 
second eigendirection (defined by ie 1 *?) is orthogonal to the direc- 
tion of the shear and the associated eigenvalue is given by '—\y\.' 
In the literature, some authors choose the opposite sign convention 
for the tidal effect term in equations and which causes an 
ambiguity in the definition for the direction of the shear between 
the two eigendirections. Our choice follows the usual weak lensing 
convention such that the direction of the image stretching distortion 
is along the direction of the shear, but implies that the direction to 
the mass that causes the tides is perpendicular to the direction of 
the tidal shear caused by it. 

The Jacobian determinant of equation Q, or equivalently, 



z - - - yz. 

z 



(2) 



This is usually r eferred to a s the Chang-Refsdal lens, after 
IChang & Refsdal <1979l Il984h . Here, y = (yi,y 2 ) T and x = 
(xi,x 2 ) T are the vectors representing the angular positions of the 
source (in the absence of the lensing) and of the image, respec- 
tively, whereas £ = yi + i>'2 and z = X\ + ix 2 are their complexified 
variables. Here, we use the overbar notation to denote the complex 
conjugation, e.g. z = X\ — ixz. In the lens equation, the effect of the 
tidal field is represented by a symmetric traceless tensor, or equiv- 
alently the complex number 7 = 7r + iyi, usually referred to as 
an external shear. The components of the (external) shear are re- 
lated to the (external) potential through 2y R = (9^ - d 2 2 )if/ ext and 

The complex lens equation can also be directly derived 
from the lensing potential by recognizing that the 'complex scatter- 
ing function' or the 'deflection function' is twice the 'Wirtinger 
derivative' (e.g. Schramm & Kavser 1995) of the lensing poten- 
tial with respect to z, i.e., the lens equation being f = z - a 
where a = 29 ; i/r. It is easy to show that this reduces to equation 
{2j with the lensing potential of the Chang-Refsdal lens given by 
if/ = (1/2) lnz + (l/4)7z 2 + C.C., where C.C. indicates the complex 
conjugate of the terms in front. In general, the further (Wirtinger) 
derivatives of a lead to the convergence d z a and the total shear d=a 
of the system. Of a particular interest is that, for a null convergent 
region (9 ; a = 0), the complex conjugate of the deflection function 
a becomes analytic I An 2005). The introduction of point masses to 
the system only adds isolated poles, and thus, the (complex conju- 
gate of the) deflection function is described by a complex meromor- 
phic (i.e., analytic everywhere except at isolated poles) function for 
any lens system without a continuous mass distribution. As for the 
Chang-Refsdal lens, we find that 



f = Z - s(z) ; 



s(z) = - + yz, 

z 



(3) 



where s(z) is not only a meromorphic function but actually a ratio- 
nal function of z (of degree 2, provided that y t 0). This implies 
that the problem of 'solving' the lens equation essentially reduces 
to finding z eros of a p olynomial (see section [3}. In addition, we 
have 9;Q- = d z a = d z s = s'(z) = —zT 2 + 7, that is to say, the total 
shear of the Chang-Refsdal lens is simply the sum of the 'internal 
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defines the area distortion factor under the linearized lens mapping. 
Hence, the reciprocal of the Jacobian determinant is (the point- 
source limit of) the ratio of the apparent solid angle covering the 
lensed image to that of the original source. Since lensing neither 
creates nor destroys photons (i.e., the surface brightness is con- 
served), this also provides us with the ratio of the apparent flux of 
the image to that of the source, in the limit of a point source. When 
the individual images corresponding to the given source position 
are not resolved, the total magnification (i.e., the flux amplification 
factor) is simply given by the sum of individual magnifications over 
all images. 

The point at which J? = is known as the crit- 
ical point, and the loci of the critical points are re ferred 
to as the critical curves (e.g. ISchneider. Ehlers. & Falcd Il992t 
IPetters. Levine. & Wambsgansdl200lh . In addition, the source po- 
sition that allows an image at a critical point is referred to as the 
caustic point, while the curves that are the mapping of the critical 
curves under the lens equation are the caustics - that is, a caustic 
point is the point on the caustics. The magnification of the image 
at the critical point and consequently the source on the caustics is 
formally infinite. The caustics are also the boundaries between the 
source positions with different number of images. (The converse is 
not necessarily true.) 

If there exist additional non-critical image positions that also 
map onto the caustics under the lens mapping, they are sometimes 
referred to as pre-caustic points. The pre-caustics 1 are the loci of 
those non-critical points th at map onto the caustics under the lens 
equation (e.g. Rhie 200*2l: lFinch et all2002l) . The pre-caustics form 
the boundaries of the image positions that are the (pre-)images of 
the interior of the caustics, whereas the critical curves are the di- 
vision lines within them between disjoint domains that the image 
cannot 'cross over'. 



Also called the transition loci by some authors 
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3 THE IMAGING POLYNOMIAL 

To find the image positions corresponding to a given source posi- 
tion f , equation J2j n eeds to be inver ted for z. The most straight- 
forward route (c.f, Witt & Mao 1995, for the binary lens equation) 
involves eliminating z from equation {2} by means of its own com- 
plex conjugate 



f(z) = l+s(z) 



■ C + - + yz ■ 

z 



1 + f z + yz 2 



This results in a rational equation of z 



g(z) = ■ 

where 

a = 

a\ = 
a 2 = 
a 3 = 



1 



■y/«-f = 



z(l + fz + yz 2 ) 



(5) 



(6) 



?y + 2\y\ 2 + |f | 2 , 
y( + 2(\y\ 2 -(, 
KM 1 - D, 



(7) 



that needs to be solved for the image positions. In general, the poly- 
nomials in the numerator and the denominator of g(z) as expressed 
in equation |6j are relative prime, and so the problem reduces to 
solving a quartic polynomial equation g n (z) = Zit=o a kZ k = 0, which 
in principle can be solved algebraically. The result further implies 
that the number of the images allowed by equation J2j for a given 
source position cannot be greater than four. 

We note however that g(zo) = is only the necessary, but not 
the sufficient, condition for zo being one of the image positions. In 
other words, not all of the zeros of g(z) are the solution of equation 

For instance, if we replace z with an unrelated new variable 
w in equations J2J and and solve for z by eliminating w, then 
the solution is still found from g(z) = for which w = f(z) is not 
necessarily the complex conjugate of z 2 Therefore, to restrict the 
zeros of g(z) to actual image positions, an additional non-algebraic 
constraint that zo = f(zo) where g(zo) = is also required. 

We also note that some properties of the Chang-Refsdal lens 
can be examined through the analysis of the rational function g(z) 
and the quartic polynomial equation g n (z) = 0. For example, 



1 - 
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(8) 



so that, if the solution zo of g(z) = corresponds to an actual im- 
age position [i.e., f(zo) = Zo], we find that g'(zo) = ^(zo). (Here 
and throughout the paper, the use of the primed symbols will be 
reserved for the total derivative.) Since the magnification must be 
real, this also implies that g'(zo) 6 R is the necessary condition for 
the zero zo of g(z) to correspond to an actual image position. In ad- 
dition, if zo corresponds to a critical point, then g(zo) = g'(zo) = 0, 
so that zo is a degenerate zero of g(z). In other words, if ( is a 
caustic point, then the imaging equation g(z) = allows a de- 
generate solution that corresponds to a critical point. Furthermore, 
g(zo) = g'(zo) = if and only if g n (z ) = g^(zo) = 0, and therefore, 
the discriminant of the polynomial g n (z) also provides an algebraic 
equation D(( c ,( c ) = for the caustics points f c to satisfy. 

We further note that, if zo is the critical point that maps to 



a cusp point on the caustics, g(zo) = g'(zo) = g"(zo) = and 
gn(zo) = g[,(zo) = g^'(zo) = 0, and thus that zo is a doubly- 
degenerate solution (i.e., a triple root) of g(z) = g n (z) = 0. In fact, 
if the deflection function is described by an analytic function, we 
can establish a connection between multiply-degenerate solutions 
of a certain analytic function derived from the deflection function 
and higher-order critical behaviour (or 'catastrophe') of the lens 
mapping at the corresponding point. However, the present analytic 
function g(z) can at most have a triple root for any ( e C so that 
the cusp catastrophe is the critical behaviour of the highest possible 
order for the Chang-Refsdal lens, provided that y 4- and |y| t 1. 

Next, for any rational function r(z) = r n {z)lr A {z) of z, where 
r„(z) and r d (z) are polynomials with no common factor, we find that 
there exist linear combinations of linearly-independent polynomi- 
als p k (z) = [rAz)f[r A (z)f- k (k = 0,...,4) that are divisible by 
g„(z). That is, Zt=o CkPk(z) = q(z)g n (z) or dividing by [r d (z)] 4 , 



E 



Cifc- 



K(z)f[r d (z)] 4 



[>d(z)] 4 



q(z)gn(z) 
Mz)] 4 



(9) 



where q(z) is a polynomial of z. Thus, for any rational function 
r(z), its values v = r(zo) at the image positions z are zeros of a 
quartic polynomial c(y) = zZt=o c * v *> provided that r(z) does not 
diverge at the image position [i.e., there is no common factor be- 
tween g n (z) and r d (z)]. In particular, if the rational function is given 
such that r(z) = R[z, /(z)] where R[z, w] is a rational function of z 
and w that satisfies that R[z, w] 6 R for w = z, then only the real 
zeros of c(v) correspond to its values at image positions. Further- 
more, if R[iv, z] = R[z, w], we can show that, for any non-real zero 
V c = r(z c ) of c(v) where Zq is a zero of g(z) that does not correspond 
to any actual image, its complex conjugate v c is also a zero of c(v). 
Hence, for this case, c(v) is completely factored into linear or irre- 
ducible quadratic polynomials with real coefficients, and therefore 
(after the common scale factor has been taken out) c(v) = reduces 
to a real-coefficient quartic polynomial equation, and only its real 
solutions are valid, which correspond to the actual image positions. 

For example, if R[z, w] = zw, then R[z, z] = zz = |z| 2 6 R 
and R[w,z] = wz = zw - R[z,w], and so, with r(z) = R[z,f(z)] = 
z/(z) = 1 + Xz + yz 2 , there exists a real quartic polynomial equation 
Zit=o t>kV k = whose real solutions are |zo| 2 where zo are the image 
positions. By equating the remainder of the polynomial division of 
2it =0 t>k(i +tz + yz 2 )* by g n (z) to be identically nil, we find the real 
coefficients (up to a common multiplicative scale factor): 
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(10) 



The result implies that, for the 4-image Chang-Refsdal lens system, 
there exists a relation among the image positions that 



1 



1 - |y| 2 



(11) 



where zi, Zi, Z3, and z\ are the complex locations of four images. 
In fact, this relation can be derived from the stricter relation found 
from the coefficients of equation Q 



2 However, if wq = /(zo) + zo where zo is a zero of g(z), then it can be 
proven that ivo is another zero of g(z) distinct from the initial zero zo. That 
is, any zero of g(z) that does not correspond to an actual image position 
always occurs in a pair. 



n 



Zk 



y 



y(\y\ 2 - 1) 



(12) 



We note that there also exist similar real coefficient quartic 
polynomial equations for 25l[zo], 23[zo], cot(argzo), c/(zo) and 
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so on Cc.f.. lAsadj2002l : lAsada. Kasai. & Kasal2 0041. In principle, 
they can be found by entirely algebraic operations outlined above 
with the choice of the rational function being z + f(z), i[/(z) - z], 
i[z + f(z)][z - /(z)]~\ or equation respectively although the 
actual algebra involved is rather messy. 

3.1 Position and Magnification Moments 

The fact that g(zo) = and g'(zo) = ^(zo) for image positions 
Zo implies that certain contour integrals involving [g(z)]~' reduce 
to the position-magni fication moment sum over th e all allowed im- 
ages JPalal & Rabmll200lt iHunter & E vans! 1200 ll) . provided that 
all zeros of g(z) are the actual image positions. Here, we show that 
this also implies that the series of the position-magnification mo- 
ment sums are coefficients for particular Taylor-series expansion of 
L?(z)] 1 (see also Evans & Hunter 2002). 

First, since lim-^ [g(z)]~' = assuming \y\ j= 1, we note that 
[g(z)] has the convergent Taylor-series expansion at z = oo such 
that 



1 _ C„ 

where h(z) = z~ 

(f) — 

JdSo g(z) 



_ l_ d" 
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1 



n\ dz" (g o h)(z) 



(13) 



z=0 



. Next, let us consider the contour integral 



dzz" 



7-i 



z"- J dz 



27tiC„ + i 




forn = 0, 1,... 
for n = -1,-2, . . . 



(14) 



where the integration path dS „ is taken to be a circle with a ra- 
dius of infinity. However, the contour integral can also be evaluated 
using Cauchy's Residue Theorem. Since g(0) t 0, non-vanishing 
residue contributions to the integral are at z = Zo where zo is a zero 
of g(z) [i.e., g(z ) = 0] 



Res 



g(z) 



= lim 



(z-z )z" 
g(z) 



?'(zo) 



(15) 



for any integer n [also assumed is that none of the zeros of g(z) is 
degenerate so that g'(zo) t when g(zo) = 0], and at z = 



Res 

z=0 



1 



z"'g(z) 



1 



1 



(m-iy. dz" 1 - 1 g(z) 



(16) 



for a negative integer n = —m (m = 1,2,...). Then, the residue 
theorem, together with equation 1 141 indicates that 



E 



g'(z k ) 



-Res z=0 {[z-"g(z)Y 



for n 
for n 



0,1,... 
-1,-2,. 



(17) 



where Zk t 0, g(zt) = 0, and g'(zt) ^ 0. Here, if all of the za-'s are 
the solutions of equation J2j - i.e., an actual image exists at every 
z = Zk - then from g'(zo) = ^fizo), the left-hand side of equation 
fPH is the same as the (signed-)magnification-weighted moment 
sum over the all images. Furthermore, the formal expressions given 
in equations 1131 and 1161 indicate that the sums in equation J 171 
are also the coefficients for the Taylor-series expansion of [g(z)]~', 
that is to say, 



— =y 

at z = <x>, or 



1 



dw m g(w- 1 ) 



. m! 7 " +1 n>( 7 



g'iZk) 



1 _ y, d" 1 
^) "Zj d^g(^) 



^ z 



?'(Zk) 



(18) 



(19) 



at z = 0. We not e that this res ult is reminiscent of Ramanujan's 
Master Theorem JBerndtiri985h . which relates moment integrals 
and Taylor coefficients. 

From the series expansion of equation JSJ at z = and z = oo, 
we can derived the following invariant position-moment sums for 
the Chang-Refsdal lens: 
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z l k f(zk) 



In addition, we also find that 
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(zk) (i - Irl 2 ) 2 ' 
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ZlJiZk) 



which subsequently implies that 



E 



j*Ci - Irl 2 ) 2 + \ 
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f(Zk) 



(20) 
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(22) 



(23) 



Similar results have als o been obtained by Witt & Mad fcOOfj) and 
IHunter & Evans! feOOll) . 



4 SOURCE ON THE SYMMETRY AXIS 

Although it is in principle possible to solve equation |6j for the im- 
age positio ns with an algebraic te chnique such as Ferrari's method 
(e.g. lAbramowitz & Stegunll972h for an arbitrary source position, 
the actual algebra is rather complicated, if not entirely uninterest- 
ing. In fact, for most cases, applying a simple numerical algorithm 
to solve the given quartic polynomial is more advantageous. How- 
ever, if the source lies along either of the symmetry axes of the 
system, the complete analytic solutions are reasonably simple to 
derive and can be illuminating. For the simplicity of the algebra, 
throughout this section, we use the lens equation J5J with the real 
direction given by the direction of the shear (i.e., the tide-causing 
mass lies along the pure imaginary axis) so that y is a positive real 
number (i.e., y = y = |y|). 



4.1 The Real Axis 

First, let us think of the case when the source is along the real axis 
(i.e., X = £). From the symmetry of the problem, we expect that 
there will be images along the real axis, too. Since any solution 
along the real axis satisfies z = z and equation {2}, we can derive 
the quadratic equation for z 



g R (z) = (1-7^-^-1=0. 



(24) 



It is a straightforward exercise to show that gR(z) is a factor of g n (z) 
if <T = X and y = y. In particular, we find that 

gu(z) = -g R (z) [r( i + r)z 2 + ( i + y)(z + y] (25) 

and therefore, the solutions of the quartic equation g n (z) = are 
r -i n 1/2 

( ±[<r 2 + 4(i- 7 )] 



Zoff 



2(1 - r) 
f , 1 / 4y 2 



2y 2y \ 1 + y 



(26) 
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Table 1. Summary of the image properties for the cases when the source lies on the real axis (parallel to the direction of the shear). Here, A = 4y - (1 + y)( 
andfi = ( 2 +4(1 -y). 



Shear 



y > 1 
y > 1 
y > 1 



Source Position 



Images 



Total Magnification M to 
71 



< y < 1 |f| < 2y(l + y)- 1 ' 2 
(K 7 < 1 |£| > 2y(l +yT l/2 
\(\<2jy~=T 



4 images [2 on (+, +) & 2 off (-, -)] 
2 images [2 on (+, -)] 
2 images [2 off (-, -)] 



4y-(l+y)^ 

(1 - y 2 )A 
m 2 - 4y + 2) 

(i-r)(-A)VS 

2y 



(y+m 

2yfy^l< \(\ < 2y(y + \y 1 ' 2 4 images [2 on (+,-)& 2 off (-, -)] l^( 4 ") / ~ 2 ~£) + . 



\& > 2y(y + l)-" 2 



2 images [2 on (-, -)] 



(y-l)AVB 

f 2 -2y 



(y-i)(-A) 



Here, two on-axis solutions z on are image locations if B = f 2 + 4(1 - 
y) > 0. Since (eR and y > 0, this happens (i) if < y < 1 or 
(ii) if y > 1 and |£| > 2 Jy—1. On the other hand, the two off-axis 
solutions z ff become physical image positions if |f | < 2y(l +y)~ 1/2 . 
The corresponding magnifications of the images can be found from 
equation J4}. We remark that the algebra may be eased by reducing 
the orders of the power of z through the use of quadratic equa- 
tions I24i and <25> for image positions and by deriving additional 
quadratic equations for magnification via a similar route using the 
resultant polynomials, as in section|3| Furthermore, since equation 
is invariant under complex conjugation, and the two off-axis 
solutions are mutual complex conjugates, they must have identi- 
cal magnifications. We find that the magnification for each off-axis 
image is 



(l+y)A 



(27) 



where A = 4y 2 - (1 + y)f 2 . If |f| < 2y(l +yT 1/2 , this is negative so 
that both off-axis images are of negative parity. If |f| = 2y(l+yT l/2 , 
the magnification diverges, implying that ( = +2y(l + y)~ 1/2 are 
the caustic points along the real axis. The resulting expressions for 
the magnifications of the on-axis images are rather complicated, 
but their sum, product, and difference, which can be derived using 
resultant quadratic polynomials, are simpler to write down; 



Mi + Mi 



\Mi - M 2 \ ■ 



?-2y 
"(l-y)A 

|l-y||A|Vs ' 



MiM 2 



1 



(l-y)AS 



(28) 



where B = f 2 + 4(1 - y) as defined following equation I26i . With 
B > 0, the two on-axis images are of the same parity if ( 1 - y)A > 0, 
whereas they are of opposite parity if (1 - y)A < 0. In particular, 
if < y < 1 and |f| < 2y(l + y)~ 1/2 , both are positive parity 
images, whilst both are negative parity images if y > 1 and |f| > 
2y(l+yr" 2 . Again, iff = ±2 yjy - 1 with y > 1 or if ( = ±2y(l + 
y)~ I/2 , the magnification of either of the on-axis images diverges, so 
that the source location corresponds to the caustic point. The total 
magnification can be found from the sum of the absolute values 
of the individual magnifications corresponding to the actual image 
positions. The results for all these case are summarized in TableQ 
We find that two negative off-axis images and one positive on-axis 
image merge into one negative on-axis image at f = ±2y(y + 1)~ 1/2 , 
whereas two opposite parity on-axis images appear or disappear 
simply by merging with each other at ( = ±2 -Jy — 1. Note that this 
reflects that the formers are actually cusp points while the latters 
are simple first-order caustic points. 



4.2 The Imaginary Axis 

The lensing of the source along the imaginary axis (i.e., I = -() can 
be analyzed in a similar manner. First, by considering the solutions 
on the imaginary axis (z = -z), we can easily factor g n (z) into two 
quadratics if \ = —£ and y = y, and subsequently find the zeros of 

g„(z): 



3[<r] + [|<r| 2 + 4(l+y)] 



1/2 



i_3[£] 

2y 



2(1 +7) 
J_(4y 2 

2y 



y 



1/2 



(29) 



Here, £ = with 3[f] e R and ( 2 = -\(\ 2 < 0. With y > 0, 

the two on-axis zeros z on are always image positions. On the other 
hand, the two off-axis zeros z off become image positions if < y < 
1 and |f | < 2y(l - y)~ 1/2 . The two off-axis solutions are related to 
each other by z <-> -z, which is also a transformation that leaves 
equation invariant, and so they have identical magnifications 
too, namely 



y 



(l-y)C 



(30) 



where C = 4y 2 - (1 - y)|f| 2 . We find M off to be positive if |f| < 
2y(l - y)~' /2 with < y < 1. The sum, product, and difference of 
the two on-axis images are similarly found to be 



Mi + M 2 = - 



Ifl 2 + 2y , 

(i+r)C 



MiM 2 = 



1 



(1 +y)CD ' 



\Mi - M 2 \ 



IgKjgP + 4y + 2) 
(l+y)|C|V5 ' 



(3D 



where D = |f| 2 + 4(1 + y). This indicates that the two on-axis im- 
ages are of opposite parity if C < and that they are both negative 
parity images if C > 0. Along the imaginary axis, the caustic points 
occur at f = +2iy(l - y)~ 1/2 with < y < 1, for which two pos- 
itive off-axis images merges with one negative on-axis image into 
one positive on-axis image. We provide a summary of these results 
including the total magnification in Table|2| 

We note that the case of the source on the imaginary axis with 
a positive real shear is actually mathematically identical to the case 
of the source along the real axis with a negative real shear, modulo 
a 90° rotation of the system. In fact, equations 12913 1 i can be ob- 
tained from equations 1261281 through the transformation y —> -y 
(and z — > iz). Nonetheless, we chose to present the case when the 
source lies on the symmetry axis that is perpendicular to the shear 
separately from that on the axis that is parallel, in order to clarify 
the physical difference between the two cases. 
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Table 2. Summary of the image properties for the cases when the source lies on the imaginary axis (perpendicular to the direction of the shear or parallel to 
the direction of the tide-causing mass). Here, C = 4y 2 - (1 - y)|<T| 2 = 4y 2 + (y - l)|f| 2 and D = \(\ 2 + 4(1 + y). 



Shear Source Position Images Total Magnification M t0 



Os?y<l |f| < 2y(l -y)- [/2 4 images [2 on (-,-)& 2 off (+, +)] 
0<y<l |f| > 2y(l -yr I/2 2 images [2 on (+, -)] 
y > 1 any 2 images [2 on (-, -)] 



4y + (1 - yU\ l 

(l-y 2 )C 
|^|(|tf +4y + 2) 

(l+yX-C)VD 
\(\ 2 + 2y 
(y+ DC 



IfJ 





0.4 0.6 

w 



2 3 

Id 



Figure 1. Top: The location of the global minimum magnification for the 
source on the real axis for < y < 1 (solid line). Also plotted are the 
locations of the caustic points (in fact, cusp points) along the real axis 
(dashed line) and along the imaginary axis (dotted line). Bottom: The ratio 
of the global minimum magnification M m j„ to the asymptotic magnification 
Masym = (1 - y 2 )" 1 at infinity. 



4.3 The Magnification along the Symmetry Axis 

For < y < 1, the total magnification of the source at the centre 
is M tal = [y(l — y 2 )Y' ■ As the source moves along the symmetry 
axis, the magnification increases until the source crosses the caus- 
tics where the magnification diverges. The caustic points on the real 
axis are located at ( = ±2y(l + y)~ 1/2 , while the ones on the imagi- 
nary axis are at ( = ±2iy(l -y)~ 1/2 . So, we deduce that the caustic is 
elongated along the imaginary axis (i.e., perpendicular to the shear 
direction or parallel to the direction of the tide-causing mass). The 
magnification initially decreases from infinity as the source moves 
past the caustic points along the symmetry axis. Along the imag- 
inary axis, it decreases monotonically and tends to (1 - y 2 )~' at 
infinity. Along the real axis, however, the global minimum magnifi- 
cation, which is slightly smaller than the asymptotic magnification, 
occurs at some finite distance from the centre, and the magnifica- 
tion increases back to the asymptotic value (1 - y 2 )~' at infinity. 
The location and the value of the global minimum magnification 
(as a function of y) can in principle be derived from the expression 



Figure 2. Magnification along the symmetry axis. Thick lines are for y = 
1 /2 while thin lines are for y = 2. The solid lines are along the real axis 
(parallel to the direction of the shear), and the dotted lines are along the 
imaginary axis (perpendicular to the direction of the shear). 

for the total magnification given in TableQ 
<rL=4y-l + — ; E = Vl-2y + 9y 2 

y 

l+y + E 1 1 - y + 4y 2 + E 
M . — '_ | 1 ' (32) 

m,n (1 -y)[l +3y 2 + (1 +y)E] \ 1 + 3y + E ' 

Figure Q shows \( mm \ and M m \ n l M-^ ym as functions of y where 
M asym = (1 — y 2 y [ is the asymptotic magnification. We find that 
Mm;,, > 1 and that M mm as M asym . In fact, for < y < 1, we have 
1 > Mnun/Masym > 9 VI/ 16 ~ 0.97428 where the minimum of 
M min /M asym occurs at y = 1/2. In Fig. [5] we plot the behaviour of 
the total magnification, as the source moves along the symmetry 
axis for y = 1/2. 

For y > 1, the total magnification for the source at the cen- 
tre is given by [2y(y + 1)]~' (which is actually the global mini- 
mum magnification), so that the lensed images are in fact demag- 
nified compared to the source. The magnification increases as the 
source moves along the symmetry axis. Along the imaginary axis, 
the magnification monotonically increases throughout and tends to 
the asymptotic value (y 2 - 1)~' at infinity. Along the real axis, there 
is a discontinuous jump in the magnification to infinity when the 
source crosses the caustic at ( = ±2 -Jy — 1. As the source moves 
across the caustic, the magnification drops from infinity to a lo- 
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cal minimum value and then diverges back to infinity as the source 
approaches the other caustic points ( = ±2y(y + 1)~' /2 along the 
real axis. Once the source moves out of the caustic, the magnifica- 
tion monotonically decreases from infinity and tends to the same 
asymptotic value (y 2 — 1) . Here, the local minimum value of 
the magnification and the corresponding source position within the 
caustic are actually easy to find. First, we note that there are four 
images - one positive and three negative parity - for the case con- 
sidered here. However, from the signed magnification sum invari- 
ant given in equation 1201 . we find that the total magnification is 
completely specified by the magnification M + of the positive par- 
ity image alone, that is, M m = 2M + - (1 - y 2 T l ■ Finally, from 
equation the minimum possible magnification of any positive 
parity image is unity when z 2 = y _1 . With y = y > 1, we find 
that i = ±y~ l/2 maps to f = +y~ I/2 (2y - 1) so that we in fact have 
2-Jy - 1 < |<T| < 2y(y + iy [/2 and the corresponding local mini- 
mum magnification is given by M lot = (2y 2 - l)/(y 2 - 1). In Fig. [2] 
we also plot the total magnification for the source on the symmetry 
axis for y = 2. 



5 CRITICAL CURVES AND CAUSTICS 

From equation t he param etric representation of the critical 
curve can be found (Witt 1990) by solving 



1 



y + e 



-2i<* 



1 + ye 



2i* ' 



(33) 



Here, the closed form solutions for the critical curve involve the 
square-root, but because of the presence of the branch-cut in the 
mapping of z l/2 , we need to be careful. Once the branch-cut of the 
mapping z 1/2 is assumed to lie along the negative real axis, the re- 
sulting expressions for the critical curve and the caustic should dif- 
fer depending on the magnitude of |y| in order to avoid crossing 
the branch-cut (i.e., there is no discontinuous jump in the resulting 
parametric expression). Hence, we find the parametric representa- 
tions of the critical curves 



(1+ye 2 *) 1 ^ (l + | r |e^) 
for < |y| < 1, and 



My 



Y^l+y-'e-^)" 2 (| 7 | + e - 2 *) 



1/2 



(34) 



(35) 



for |y| > 1. Here, y = |y|e 2 "^ and <j> = <p - <p r . The parametric 
representation of the caustics can be found by mapping the critical 
curve though the lens equation 



Zc(0) 



(36) 



However, before we actually derive expressions, let us exam- 
ine first 



= e* \^*z'M) + e^(0)l = 2^%[e-^z' c (<P)l 



(37) 



In other words, we find that <T t '(0) is parallel to e 1 *, that is, the pa- 
rameter (p is defined such that it is the argument angle of the tan- 
gent to the caustics at the caustic point. Furthermore, we also find 
that the caustics form a cusp point [££(0) = 0] if e~ 1( *z^(0) is pure 
imaginary number, that is z' c (<f>) = ire 1(S with rel. This indicates 
that 4(0) is orthogonal to e 1(S at the critical point corresponding to 
a cusp point. In fact, the general condition for the cusp along the 



caustics is that the tangent to the caustic is the normal direction to 
the critical curve at the corresponding caustic point and the critical 
point. 

Now, the actual parametric expression of the caustics of the 
Chang-Refsdal lens are given by 

1 1 + 2ye- 21 * 



U4>) : 



(1 +fe 2i *) 1/2 (l+ye- 21 *) 
1 1 + 2|y|e 



1/2 



/ -\l/2 / -\l/2 

(l + |y|e 2 '<*) (l + lyle- 21 * 4 ) 



(38) 



for < |y| < 1, and 



1 



1 



-lyl- 



2 + y-'e- 



l P 2ic* 



(1 +r -i e -2. ( 6) 1 / 2 (l+y-le 2 ^) 1 ' 2 
1 2|y| + e 2i * 



(|y| +e- 2i *)' /2 (|y| + e 2 *) 



1/2 



■Mr 



(39) 



for |y| > 1. Here, we note that the critical curves and the caustics are 
single simply-connected curves with <p e [0, 2n) for < |y| < 1, 
while they are two separate mirror-symmetric closed curves with 
(f> 6 [0, 7t) f or y > 1 . 

The cusp points along the caustics can be found using 



1 



CW = 2^3 

_(1 + ye- 2 *) ' 
for < |y| < 1, and 



(y + e 2i<?) 3 / 2 



(40) 



(41) 



for ]y] > 1. First, for < |y| < 1, we find that the condition C(<t>) = 
is equivalent to (1 + |y|e 2l *) 3/2 6 R so that we have 2i0 = nn for 
corresponding to the cusp points. Here, n is an integer. Hence, the 
caustics have four cusps at 

2M , ..-wv 2 lrl 



fix 



(42) 



(1 + lyl)" 2 ' d-lyl) 1/2 ' 

which is consistent with the result in the previous section. Next, for 
|y| > 1, the condition that CW = reduces to (|y]e 2 '* + l) being 
a positive real number, or equivalently arg(l + |y|e 2 '*) = (2«7t)/3 
where n is again an integer. From simple geometric considerations, 
this reduces to 20 = 2«7t, 2rm + (27t)/3 + f, or 2rm + (4n)/3 - f 
where £ = arcsin[ V3/(2|y|)]. Therefore, we find that there are three 
cusps along each caustic (hence, in total six, 2 x 3 = 6) at 

V2e'f /2 



±e'' 



Cc = ±e" 



(V4|y| 2 -3- l) 

V2e-'^ 2 
(V4|y| 2 -3- l) 



-i 2 -3 
2\y\' 12 



-(2- V4W T3 3 + e 2m/3 ) 
(2- V4W Tr 3 + e- 2m/3 ) 



1/2 



(i + lyh 1 ) 



lW 2 ' 



(43) 



5.1 The Areas under the Critical Curves and the Caustics 

Next, since e 2l(4 is periodic for <p (with period 7t), one can conclude 
that z c ((p) [and consequently £ c (0)] is also periodic although the 
period is not necessarily n. If there exists a convergent Fourier- 
series expansion 



»(0 = £ 



\k(s)<p . 



C(<f>) = J] dke 

k=-oo 



\kvjd) 



(44) 
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where oj is the angular frequency of z c ((p) [i.e., the ^-period of z c (</>) 
is In/co], then from 



(l-ku>)<p 



kc k + \k c 2/u _, 



(45) 



we find that kd k = kc k + (k—2/(o)c2/a>-k, assuming 2/lj is an integer 
[i.e, the 0-period of z c (0) is an integer multiple of n]. Here, the 
relations for k = and k = 2/a> together indicate that C2/ a = = 
0, whereas c and d are unspecified. Then the caustics may be 
written in terms of c k 's as 



fe(0 - d « = Z 



k-2/o>_ 

Ck H ; C2/<j-t 



Zj 

fcS2/w 



(46) 



One application for which the Fourier-series expansions may 
be useful is to find the area within the critical curve and the caus- 
tics. We note that the 'area 2-form' in C may be written as dhc = 
d(zdz)/(2i). 3 Then, the (signed) area within a simply-connected 
closed curve dC = [z(p)\z(p) : [0, P)^> C,pe [0, P) c R) where p 
is a parameter and P is the period can be found to be, from the 
fundamental theorem of multi variate calculus - the ' generalized 
Stake's theorem' (c.f.. lGould & Gaucherellll99l lDominiklll998h 



S = fd 2 x = ^ f d(zdz) = i £ 1 dz = ^ f 
Jc 2i J c 2i Jg C 2i J 



zip) z\p)Ap. 



(47) 



The area within the critical curve can be expressed in terms of its 
Fourier coefficients, from the generalized Parseval's theorem, as 



sum 



r 2n/w vh 

Zc(0)4(0)# = n )k\c k \ 2 
Jo *-r 



(48) 



Provided (i) that the period of f o (0) is a ^ so 2n/o> and (ii) that the 
caustic <f c (0) does not exhibit self-intersection, the area within the 
caustic can also be written in similar manner 



i Ml * \ 



■2/w-k 



(49) 



k*2/K. 



We note that this is typically negative as the orientation chosen by 
the parameter <f> is usually clockwise. 

For < |y| < 1, from equation I34i . we find the Fourier- series 
expansions of the critical curves and the caustics to be 



z,c 



...TO*"™" 



S 



(50) 



fer*r = £ y [ (2it - De (2t+I) * + (2k + l)e- (2 *- lW ] |y|* (51) 



where (a)j is Pochhammer symbol such that 

T(a + it) 



(a)^ = oia + 1) • • • (a + k — 1) 



r(a) 



(52) 



d(zdz) = dz A dz = (dx\ - idx2) A (dx\ + id;t 2 ) = 2i dxi A d*2 



and 



C k = (-!)' 



H*- 1/2) 

fe!T(l/2) ' 



(53) 



We note that both series expansions are also Taylor-series expan- 
sions with respect to y at y = 0. If we truncate the expansion of 
the caustic after (k = l)-term, the expression reduces to the equa- 
tion of the tetra-cuspi-hypocycloid, or the astroid. The fact that 
the caustic for the Chang-Refsdal lens with < \y\ < 1 has an 
astroid-like shape has been widely acknowledged. In fact, the as- 
troid is the generic shape of the causti c for a lens syst em under a 
weak quadrupole perturbation jKovneJl987t lAr]|2005h . Here, the 
0-period of both curves are 27t, and both are simply connected. 
Thus, the area under each curve is 



S[z c ] = Yj 



2T(k+ l/2)r(fc + 3/2) 2k 
T(k + if 171 



n,i\\~,\,l-,\y\ 2 



■■ 2n 



2T(k+l/2)T(k + 5/2) 2k+2 
T(k + 2) 2 171 



(54) 



(55) 



where 2 F\ (<z, b\ c\ x) is the Gaussian hypergeometric function. The 
result for the area under the caustics expressed using the elliptic 
integrals (see Appendix[Aj can be compared to the area und er the 
planetary caustics due to an 'external' planet (see e.g. lBozz^ 2000). 
which is basically in the same form up to the overall scale factor 
(oc q) once \y\ oc d~ 2 has been properly chosen, where q and d 
are the mass ratio and the separation in units of Einstein ring of 
the planetary companion. This is because the effect of the primary 
star on microlensing due to the plan et can be w ell approximated 
by the Chang-Re fsdal lens iGould & Loeblll992tlGaudi & Gould 
1997; Bozza 1999). This is also the case for any point mass lens 
perturbed by mass lying far outside of its Einstein ring radius, as in 
a wide- se paration binary for instance toominikt 1999; Bozz al^OOfl : 
lArl2005l) . 

For |y| > 1, we find the Fourier- series expansion of the critical 
curve as 



Zc<Z 



-\y\w\ n k \ \2j h |y|* j 



(56) 



and the caustic 



= +lrl 1/2 j do + £ c k (k 1)6 ,7 } ■ ' 57 - 



where C k is given by equation 1531 and 



2- ■ 



Irl' 



(58) 



Unlike the previous case, the leading terms are given by a con- 
stant. Hence, we find that the critical curves are 'centred' at z = 
±|y|-'/ 2 e i ^, and that the caustics are at f = ±(]yr 1/2 - 2|y| 1/2 )e'^. 
The truncation of the series for the caustic after the first non- 
constant term for this case leads to the equation of the tri-cuspi- 
hypocycloid, or the deltoid. In other words, as |y| — > 00, the caustics 
reduce to a pair of mirror- symmetric deltoids. Here, the ^-period 
of each separate curve is given by n so that the area under each is 
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Figure 3. Critical curves, pre-caustics, and caustics of the Chang-Refsdal 
lens. The top panels are for y = 1/2 and the bottom panels are for y = 2. 
The left panels show the critical curves (thick lines) and the pre-caustics 
(thin lines), and the right panels show the caustics. The numbers indicate 
corresponding points mapped under the lens mapping. 



found to be 



T(k + 3/2) 2 



j r(* + i)r(* + 2) | r p 



7t /3 3 J_ 

4| 7 p 2 ! U' 2 ; ; w 2 



^2 2 

r(* + 3/2) 2 



1 



fc=0 
■K 



T{k + i)r(* + 3) iri 2 



,3 3 1 



2' 2' 



(59) 



(60) 



Here, the area derived is for each one of the curves, and therefore, 
the total area under the critical curves and the caustics is twice this 
result. Again, the results expressed in terms of elliptic integrals are 
provided in Appendix lAl 



5.2 Pre-caustics 

Since the critical points are the degenerate zeros of the imaging 
polynomial corresponding to the source on a caustic point, the cor- 
responding pre-caustic points can be found from the remaining ze- 
ros of the polynomial after the two degenerate linear factors (z-z c ) 2 
have been factored out. The coefficients of the resulting quadratic 
polynomials piz 2 + piz + po are (for simplicity, y is taken to be real) 



P2 

Po 



rd-r ); 

1 + ye 2 '* 



-ye- 2 '""(l + ye 2 "* 1 ) 



(l-2y 2 )(l+ye 2i *) 



1/2 



(61) 



pi = r(r - !) 
Pi = 



i i o,vj\i/2 1 + ye- 
(2y 2 - l)(y + e- 2,(S ) + 



2i<» 



(y + e 2i *) 1/2 



y(y + e- 2 '*) 



(62) 



for y > 1. Finally, we find the parametric expression of the pre- 
caustics to be 



(l + ye 2 ^) 



1/2 



Z P 6 2y(l - y 2 ) 



(\ +ye 2i,A 1/2 + r 

1 - 2y 2 - 



1/2 



2y(l - y 2 ) 



1 - 2y ; 



(1 +ye- 2i *) 1/2 
(l+ye 2 '*) 1/2 -r 



(1 +ye- 2i <*) 1/2 



for < y < 1 , and 

(y + e- 2 *)" 2 
Zp ~ ± 2y(y 2 - 1) 

for y > 1 . Here, 



2y 2 



1 + 



(y + e- 2 ^) 



1/2 



(y + e 2i -») 1/2 



r i i i \ 1/2 

— = 1 + y cos 2<f> - ( 1 - 2y ) 1 1 + y + 2y cos 2(f> 1 



(63) 



(64) 



(65) 



(66) 



We note that, for < y < 1, the pre-caustics are two separate 
curves, one of which zt(<f>) completely encloses the critical curve 
and the other of which z„((f>) is completely enclosed by the critical 
curve. We also find that the curve z„((f>) is in the generic shape of 
the bi-cuspi-epicycloid, or the nephroid whereas zjj(0), as a whole, 
appears to be the inverse curve of z p (</0 with respect to the ori- 
gin, rotated by 90° and rescaled by a factor of -^/l - y 2 . On the 
other hand, for y > 1, the pre-caustics are two mirror-symmetric 
curves, each of which encloses one of the mirror symmetric critical 
curves. For both cases, the critical curves and the pre-caustics are 
co-tangent to one another at the critical points corresponding to the 
cusp points. The examples for y = 1/2 and y = 2 are shown in 
Fig.0 

Following Fi nch et all J2002t) . the ratio of the areas under the 
caustics and the pre-caustics provides us with the mean total mag- 
nification of the source that lies within the caustic. Although the 
complicated parametric forms of the pre-caustics indicate that the 
analytic evaluation of the area under them using the present ex- 
pression is tough, following equation J47> . it reduces to a one- 
dimensional (complex-contour) integral, which is straightforward 
to evaluate numerically. In fact, it turns out that the area under the 
pre-caustics can be evaluated analytically. The result for y > 1 is 
particularly notable and is given in the next section. On the other 
hand, the result for ^ y < 1 is more complicated and is simply 
quoted in Appendixlcl 

It is known that the asymptotic form of the magnification dis- 
tribution is given by P(M) ~ M~ 3 if the effect of the divergent 
magn ification in the vicinity of the caustics is dominant l Schneider 
1987). If the full magnification distribution for the source in the 
caustics is well approximated by its asymptotic form P(M) ~ M~ 3 
down to the minimum possible value M c . , the mean magnification 
is given by (M) = 2M c min . Here, M^ in = [y(l - y 2 )Y l for < y < 1 
and M L mm = (2y 2 - l)/(y 2 - 1) for y > 1. In Fig.H we plot <M> and 
(M)IM C . as a function of |y| to examine how much the magnifica- 
tion distribution deviates from its asymptotic form. Also plotted in 
Fig.|5|are the effective power-index n cff for the magnification dis- 
tribution defined such that (M) = [(n c ff - l)/(ieff - 2)]M^ in , which 
would be the mean magnification if the magnification distribution 
were characterized by the power-law P oc M~" dl , and the effective 
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Figure 4. Top: the mean total magnifications (M) for the source within 

the caustic, that is, when there are four images, as a function of \y\. Also 

shown as a dotted line is the minimum possible magnification M c . for the 

r & min 

same configuration (that is, 4-image system). Bottom: the ratio (M>/M c . . 
The horizontal dashed line marks the ratio of 2, which corresponds to the 
case that the magnification distribution is given by a truncated power law 
PocM' 3 . 



4th-order coefficient f eff for the magnification distribution given 
from (M) = [3(2 + c cff )/(3 + 2c e fj)]M^ in , which would be the mean 
magnification if the magnification distribution were characterized 
by P(M) oc M- 3 (l+c cff AT m /M) (note that 1.5M^ im < (M) < 3AT in 
for c e ff > — 1, that is P(M) being non-negative). We find that, for 
\y\ > 1.1, (M)/M c mm ss 2 within 5%, implying that P(M) ~ M~ 3 
may indeed be a good approximation for the magnification distri- 
bution for the source within the deltoid-like caustics down to the 
minimum possible magnification. In fact, we find the full analytic 
expression for the magnification distribution for this case in section 
|6|and show that it is well approximated by P(M) ~ M~ 3 . On the 
other hand, for < \y\ < 1, (M)/M^ m is smaller than 2, which im- 
plies that the power-law index should be steeper than -3 if the dis- 
tribution were a truncated power-law, or that there is a significant 
positive higher order term if P(M) is approximated by a power- 
series to M with the leading term as M — > oo given by oc M~ 3 . We 
also note that the limiting value of (M)/M^ ia as y — > is found 
to be 8(37t)- 1 [3K(l/2) - 2E(l/2)] « 1.80149, which translates to 
approximately c eS ss 1. Here, K(x) and E(x) are complete ellip- 
tic integrals of the first and second kind (see Appendix IaI for the 
notation). 



6 LINES OF EQUAL INDIVIDUAL MAGNIFICATION 

From equations and I33i . it is not only the critical curves but 
also the lines of equal magnification for an individual image that 
can be written parametrically. That is, the solution of 

— - y = mer n * (67) 

7 

defines the locus of image positions for which the magnification 
is given by 1 1 — m 2 \ 2 . The area under each curve is proportional to 



Figure 5. Top: the effective power-index n^g for the magnification distribu- 
tion. Bottom: the effective fourth-order coefficient c e fj in the magnification 
distribution as a function of \y\. See the text for details. 



the cumulative probability that a given image is at least magnified 
by 1 1 — m 2 | 2 , and thus, its derivative with respect to m is related 
to the magnification distribution for a given image. Although this 
is not the distribution of the total magnification for a given source 
position that we want, it is still of some interest to examine this 
distribution. 

As before, we examine different ranges of |y| separately, but 
now the dividing line occurs at \y\ = m in the general case. First, 
we consider the case that < \y\ < m. For this case, the line of 
equal magnification is 



,,11 + Me* 



k\ 



n k+l/2 ' 



,(2t+l)i0 



(68) 



This is a single simply-connected curve. The areas under it and its 
derivative are 



S[z m ] = £ 



2r(fc + 3/2)r(fc+ 1/2) \y\ 2 
T(k + l) 2 m 2k - 
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(69) 



(70) 



The derivative can simply be transported to the source domain by 
multiplying by the Jacobian determinant (i.e., the inverse magnifi- 
cation), namely 



dm 



\m 2 - 1| 



dS[z m ] 



Am 



n\m 2 - 1| 



3 3 



Irl 2 



Fl «2'2 ;1; ^ ' (71) 



The 'area' under the curve 
(2k- l)e <2 * + '« 
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2m k+l ' 2 



2m*-i/2 



lyl' 1 



(72) 
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(where C* is given by equation l53l . is in principle found either by 
integrating this or by a similar method as in section [5] the result 
however is not necessarily the actual area since the curve f m may 
self-intersect, but the signed sum of the areas of adjacent simply- 
connected regions. 

If |y| > m > 0, we find that the loci of the equally magnified 



images 



1 / m _ 
- — T7T 1 + — e 



-1/2 



(-If 

k\ 



Irl 



t+l/2 



(73) 



are two separate mirror-symmetric simply-connected curves. The 
area and its derivative for one of the simply-connected curves is 



S[z m ] = £ 



r(/t + 3/2) 2 m 
r(k + i)r(k + 2) i^F 



Tim 1 3 3 m" 
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whereas 
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(74) 



(75) 
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(77) 



Now, if |y| > 1 > m > 0, equation J76I provides us with the dis- 
tribution of the magnification of all positive parity image(s) for a 
given source position because there is at most one positive par- 
ity image for any source position for |y| > 1. More importantly, 
the total magnification accounting both positive and negative parity 
images when the source lies within the deltoid caustics is given by 



M, 



+ 



1 



Irl 2 - i ' 



\ — m2 

and thus, 

dS[<f m ] dm 
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(78) 
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Since < m < 1, we have 
2|y| 2 - 1 



(79) 



(80) 



(81) 



Irl 2 - 1 

While equation 1791 is an exact analytic expression for the con- 
ditional distribution of the magnification for a source within the 
deltoid caustics of the Chang-Refsdal lens with |y| > 1, the elliptic 
integrals (to which the hypergeometric function reduces; see Ap- 
pendix [A) are not the most convenient functions to deal with, and 
therefore, it may be better to come up with a simple and reasonable 



approximation. First, we note that the leading term of the Taylor- 
series expansion of equation 1791 at M tot = oo is found to be 



dS[C„l 



dM, 



Irl 3 



3 3 
2' 2' 
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IrP 
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(82) 



which confirms the general result on the asymptotic form of the 
magnification distribution (~ M~ 3 ). Next, we consider the limit- 
ing value at M c mm = (2|y| 2 - l)/(|y| 2 - 1), that is, d5 [f,„]/dM lot = 
7t/(2|y|) 3 . We find that the coefficient for the leading term of the 
Taylor series is approximately equal to 7t[M^ in /(2|y|)] 3 for y > 1.1. 
In other words, the approximation of the distribution by a simple 
truncated power-law 

P Wot )«P(M^)(^l) ®(M tot -M^J (83) 

V min / 

extending down to the minimum magnification (eq. 1811 is in fact a 
reasonably good approximation for y > 1.1. 

Since we have the expression for the full magnification dis- 
tribution, the mean magnification for the source within the deltoid 
caustics may be derived from the normalized moment. However, 
the same special property of having one positive parity image to- 
gether with the existence of the magnification invariant indicates 
that it can be derived from the ratio of the areas under the caustics 
and the critical curve. That is to say, for sources within the deltoid 
caustics, their corresponding positive parity images completely fill 
the region inside the critical curve, and thus, the ratio of the area 
under the critical curve to that under the caustic 

2 2^(3/2, 3/2; 2; | 7 |- 



<M + > = 



2 ) 



(84) 



2^(3/2, 3/2; 3; |y|- 2 ) 

gives the mean magnification of the positive parity image. Finally, 
equation 1201 implies that the mean total magnification (M tot ) for 
this case is related to the mean magnification of the positive parity 
image (M + ) through 

4|y| 2 2^1(1/23/2; 2; |y|- 2 ) 



<M tot > = 2<M + > + 



|y| 2 -l |y| 2 - 12^(3/2,3/2; 3; |yh 2 )' 



(85) 



Using this result, we can also derive the analytic expression for the 
area under the pre-caustics in equation <65t as 



S[z p ] = {M tot )S[( c ] 



2|yl(lyl 2 -D 



±^;2; J- 
2 2' '|y| 2 



(86) 



7 A CONVERGENT BACKGROUND 

One possible generalization of the Chang-Refsdal lens equation J2j 
is the addition of a constant convergence term 

1 



■ JZ - KZ, 



(87) 



which describes a continuous mass distribution with a constant sur- 
face density. While the effect of the constant convergence term (kz) 
is basically an overall magnification (uniform defocusing) of the 
lens domain with respect to the source domain, the specific be- 
haviour of the lensed images is quite distinct, depending on the sign 
of (1 - k). If < k < 1 (note that k must be non-negative and real 
since Poisson's equation indicates that k is simply the rescaled sur- 
face mass density), introductions of new variables rj = (/(I - k) [/2 , 
£ = {I - k) 1/2 z, and g = y/( l - k) simply red uce the lens equation 
to the form of equation (Paczvnski 1986). Hence, the lensing 
behaviour of these systems basically follows that of the standard 
Chang-Refsdal lens discussed so far, except for the use of the 're- 
duced shear' g = y/(l —k) and the additional focusing factor (1 - k) 
in the magnification. 
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Figure 6. Caustics for the Chang-Refsdal lens with additional background 
convergence in the cases g = 1/2, V3/2, 0.95, and 2. 



However, for k > 1, which is sometimes referred to as 
the over-focusing case, the introduction of the rescaled variables 
leaves the sign of the point-mass deflection term changed. With 



r, = -£/(*- l) 1 ' 2 , 
lens equation 

1 = £ + | + sl 
At the limit of g = 



g = (k - 1) ' z, and g = j/(k - 1), we have the 



s 2 



l l'/l 2 



0, the system allows two images 

1/21 



(88) 



(89) 



if \t]\ > 2, whereas no image is formed if \rj\ < 2. The total magnifi- 
cation for > 2 is given by 



(1 - k) 1 M k 



\n\ 2 



(90) 



where the (1 - k) 2 factor is due to the overall focusing factor. Note 
that the magnification diverges at \rf\ = 2 indicating that the caustic 
of this system is given by the circle \rf\ = 2. The corresponding 
critical curve is found to be another circle |f | = 1. For this system, 
when the source crosses the caustic into its interior, two images, 
one inside and the other outside of the critical curve divergently 
merge and vanish. Two images re-appear when the source crosses 
the caustics outwards. 

For the general case when g # 0, the lensing behaviour can be 
studied using the same methods developed throughout this paper. 
For example, the image positions can be found by solving similar 
rational equations or polynomials as equations J6j and 0, but with 
the following replacement of variables: z — > f, y — > g, y — > g, 
f — > 77, and I — > —fj. One notable difference between the new 
polynomial equation as compared to that of the standard Chang- 
Refsdal lens is that it has a quadruple zero if g = VI/2 and 277 = 
±(3 + VS) ± (3 - V3)i, which indicates that this system will exhibit 
third-order critical behaviour, namely the swallow-tail catastrophe 
along the caustics, for g = y3/2. In fact, this result can be more 



easily obtained by analysis of the lens mapping along the lines of 
section[5] 

In Fig.|6| we plot the caustics for these systems with g = 1/2, 
V3/2, 0.95, and 2. The topology of the critical curves are basically 
the same as the standard null-convergent Chang-Refsdal lens. If 
< g < V3/2, the caustics are the oval-shape cuspless curve, 
within which no image is formed. At g = v3/2, the caustic be- 
comes a pin-cushion shape with a swallow-tail catastrophe at its 
four corners. Again, the system allows two image outside the caus- 
tics and no image inside the caustic. For V3/2 < g < 1, the 
swallow-tail catastrophe of the caustics metamorphoses to self- 
intersecting curves with two cusps at each of four corners. There 
exist four images for the source within the fish-tail-like region of 
the caustics, while there is still no image for the source within the 
main pillow-like region of the caustics. As for g > 1, there appear 
two deltoid-like caustics, but unlike the standard null-convergent 
case, the caustics are stretched in the length direction rather than 
the width direction. 

While a more detailed study of the lensing behaviour can be 
easily done, it is beyond our scope here. However, we note that such 
a study may help understand the distinct microlensing behaviour of 
the central demagnified macro-image at the maximum of the time- 
delay surface (see e.g. Dobler, Keeton, & Wambseanss 2005). 



8 N POINT MASSES WITH SHEAR 

Another physically interesting generalization of the Chang-Refsdal 
lens is to the case of N point masses with shear. This can be repre- 
sented by the lens equation 



V m ' 



■ yz, 



(91) 



where m, are the relative masses and /, are the complex positions. 
We find that the deflection function (c.f, eq.|3J 



Z - J z - h 



yz : 



yznfiz-id + ^m.^iz-lj) 

nffe-w 



(92) 



is a rational function of degree (N + 1), where N is the number 
of point masses. To establish the imaging equation, we follow the 
same procedure as in section|3] This leads to consideration of the 
rational equation g(z) = z - h(z), or equivalently the fixed points of 
the mapping h(z) = (f°f)(z) where f(z) = £+s(z)andf(w) = f(w). 
Here, since f(z) is of degree (N + 1), h{z) is of degree (N + l) 2 . 
In addition, as z — » 00, s(z) [and consequently /(z)] diverges lin- 
early (7 t 0) and so does h(z). This indicates that the degree of 
the rational function g(z), as well as the order of the polynomial 
in its numerator, is (N + 1)- <Wittlll990j) . Obviously, as N grows, 
finding image positions through solving (N + l) 2 -th order polyno- 
mial equation quickly becomes impractical. On the other hand, the 
parametric representation of the critical curves (and subsequently 
the caustics) can be derived from |s'(z)| 2 = 1 or by simply solv- 
ing s'(z) = -e~ 2l< *, which reduces to a (2/V)-th order polynomial 
equation. Here, the parameter <p is chosen such that cf> is the argu- 
ment angle of the tangent to the caustics (section[5}. In other words, 
finding the critical curves and the caustics is comparatively easier 
as the order of the relevant polynomial grows linearly with N. By 
contrast, solving the lens equation requires us to solve a polynomial 
whose order grows quadratically as N increases. 

One additional point to note regarding the N point mass lenses 
under constant external shear is that, despite the order of the imag- 
ing polynomial being (N + l) 2 , the upper bound for the number 
of images grows linearly with N. Using the result of Appendix 
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IbI we find that the total number of images for a lens system 
with N point masses with non-zero shear may not be greater than 
5(N + 1) - 6 = 5N - 1. [Compare th is to the upper bo und for 
7 = 0, which is 5(N - 1) with N > 2 <Rhidl200lU2003l) and for 
which the order of the corresponding imaging polynomial is N 2 + 1 
.] This also means that the imaging polynomial must 
have a zero that is not a image position if N > 3 (for y t and 
|y| + 1). Furthermore, if < |y| < 1, we find that n + < 2N and 
«_ — n+ = N - 1 where n + and n_ are the numbers of positive 
and negative parity images, respectively. For \y\ > 1, these results 
become n + < 2N — 1 and »_ - n+ = N + 1. 

We note that, for the least extreme case (i.e., the source 'out- 
side' the caustics), the difference in the number of opposite parity 
images can be understood as being due to one negative parity image 
associated with each point mass, plus one additional image of ei- 
ther positive (0 < |y| < 1) or negative (|y| > 1) parity. Here, the one 
additional image of 'spare' parity can be regarded as the weakly 
perturbed macro image: the system with < |y| < 1 represents the 
macro image at the 'minimum' in the time-delay surface (the posi- 
tive parity image at sub-unity convergence), while that with |y| > 1 
corresponds to the macro image at the 'saddle point' in the time- 
delay surface (the negative parity image). 



9 CONCLUSIONS 

This paper has presented a complete mathematical treatment of 
the Chang-Refsdal lens. This simple and flexible model has found 
widespread applications in astrophysics. The lens equation for the 
Chang-Refsdal lens can be converted to a polynomial of the fourth 
degree, the roots of which provide the image locations (together 
with possibly spurious roots). If the source lies on one of the sym- 
metry axes with respect to the external shear, this quartic readily 
factorizes into a product of quadratics, which enables the lens equa- 
tion to be solved via elementary means for these special cases. 
Simple lensing invariants exist for the products of the (complex) 
positions of the four images, as well as for sums of the moments 
of their signed magnifications. In other words, provided the source 
lies within the astroid or deltoid caustics, the sums of signed mag- 
nification, multiplied by moments of the complex positions, are in- 
variant. 

Of the results in this paper, we highlight just two here. First, 
we have provided the equations for the pre-caustics. Only one 
other lens model has been known for which the pre-caustics can be 
writte n down, nam ely the singular isothermal sphere with external 
shear iFinch et aljEool . Second, the exact analytical expression 
(eci. l79> in terms of hypergeometric functions (or elliptic integrals) 
has been calculated for the distribution of magnifications for the 
source within the deltoid caustics of the Chang-Refsdal lens. For 
many practical purposes, however, we have shown that a truncated 
power-law (eq. l83> is an excellent approximation to this distribu- 
tion. There are a number of applications of these two results. For 
example, the Chang-Refsdal lens is often used as an approximation 
to the binary lens equation, especially in the case when the lens is a 
star plus planet. The areas under the critical curves and pre-caustics 
are directly related to the mean magnifications of the 2 and 4 im- 
age geometries, and hence the detectability of a planet through the 
gravitational microlensing effect. Equally, the distribution of mag- 
nifications is applicable to the statistics of high-magnification mi- 
crolensing events in the lightcurves of quasars. 
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APPENDIX A: ELLIPTIC INTEGRALS FROM 
HYPERGEOMETRIC FUNCTIONS 

We note that all hypergeometric functions found in the main text re- 
duce to expressions involving elliptic integrals. Although the gen- 
erality of the hypergeometric function can be seen as grounds for 
superiority, there are some advantages in using elliptic integrals es- 
pecially for computational purpose, partly because numerical rou- 
tines to evaluate elliptic integrals are more readily available than 
hypergeometric functions (see e.g. lPress et alll992li . For this rea- 
son, and also to help readers compare results to those found in the 
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literature, we provide alternative expressions written in terms of 
elliptic integrals. 

First, we note a relationship between hypergeometric func- 
tions with parameters differing by unity; 



2F1 {a, b; c; x) — jF\ {a - 



l,b; c; x) = —X2F1 (a,b + 1; c + l;x), 
c 



and one of the Euler transformations of the hypergeometric func- 
tions; 



2F1 (a, b; c; x) = (1 



x) c -"- b 1 F l (c -a,c- b; c; x) . 



Using these two, we can express any Gaussian hypergeometric 
function 2F1 (a, b; c; x) with a and b being half-integers and c a pos- 
itive integer in terms of a linear combination of 2 Fi(l/2, 1/2; l;x) 
and 2F1 (—1/2, 1/2; l;x) with the coefficients being rational func- 
tions of x. However, the complete elliptic integrals of the first and 
the second kind are actually hypergeometric functions 



K(k) 



E(k) 



d<p 



^ ^1 - k 2 sin 2 

f ' d<pjl^ 

Jo 



k 2 sin 2 



nil 

2'2 :1; * 



7t Ml , 
2 2F 'h'2 ;1 ^ 



In other words, the hypergeometric function 2 F\ (a, b; c; k 2 ) is in 
general expressible in terms of the complete elliptic integrals if a 
and b are half-integers and c is a positive integer. For example, 



I I 

>2'2 



; i;k' 



2 \ 22 



1 1 

~2'~2 



; l;k' 



W) 
l-k 2 

2E(k) - K(k) 
2E(k)- (1 -k 2 )K(k) 



etc. Fur ther example s can be found in standard refere nces 
such as lErdebdeLalJ J^QQd), lAbramowitz & Stegunl 1 19721) . or 
Gradshtevn & RvzMy 11994 . or on-line resources like Wolfram 
Research site 4 . (Be careful as many authors use various different 
conventions for the arguments of the elliptic integrals!) 

From these, we find alternative expressions written in terms of 
elliptic integrals for equations J54I and J55I : 



S[z c ] 



2E(lrl) 



1 - \y\ 2 

S[fiJ = 27t + 4K(M)-8E(M), 
for equations 1591 and 1601 : 

5Uc] = w^t e (r) _ R k (r 

2(2|y| 2 - 1) / 1 \ / 1 

for equations |69}, OHJ, and J7TJ; 



S[z m ] = 



2m 



m 2 - |y| 2 

dS[z,„] = 2_ 

dm 



Irl 



d/77 



■** - |y| 
2|m 2 - 1| 



Kl^l- 



lyl 2 



m , 
2m 2 



Irl 



Irl 2 



m 2 - |y| 2 \ ra 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 
(A6) 
(A7) 



http://functions.wolfram.com 



for equations 1741 . 1751 . and 1761 

Irl 



S[z m ] = 

dS[z,„] 
d/;j 

dS&] 



Irl 2 



m 
m 



El— 1-— Kl — 

2 Mrli Irl I Irl 



2|yl 2 



dm 



lyKIrl 2 - « 2 ) 

m\l — m 2 \ 

: IrKIrl 2 -™ 2 ) 



|y| 2 - m 2 
2|yl 2 



m \ I m 

M/ 1 \M 



Irl 2 



— I - K| — 



Irl 



\lrl 



and for equation 1861 ; 



2lrl 



Irl 2 



1 



K — 

Irl 



E — 

Jrl 



(A8) 
(A9) 
(A10) 

(All) 



In the literature, sometimes results are given in terms of the 
incomplete elliptic integrals even if they are entirely expressible 
using only the complete elliptic integrals. Let us consider the Jacobi 
form of the incomplete elliptic integrals of the first and the second 
kind; 



F(s;k) 
E(s;k) 



r * 

Jo (1 - v 2 )'/ 2 (l - £V)'/ 2 

r d-fc 2 v 2 ) i/2 , 

J irr^iTr*'- 



Here, we note that K(/c) = F(l ; k) and E(k) = E(l; k). Then, we find 
that, with s = l/k = z' 1 . 



F.I; 



E|-;z 



dv 



(1 _ v 2)l/2 (1 _ z 2 v 2)l/2 



dw/z 



Jo (1 - u' 2 /: 2 ) 1 -(! -ii- 2 ) 1 - 

h (1 _ z 2 v 2 } l/2 



Ml 

z \ z 



T 1 (l-w 2 ) 1/2 dw 
V "Jo (l-w 2 /z 2 )'/ 2 T 



(1 - v 2 )" 2 
l-z 2 + z 2 (l-w 2 /z 2 ) dw 

o (i-w 2 )" 2 a-w 2 /z 2 y 12 z 

so that z^K(z-') = F(z~';z) and zE^" 1 ) = E( Z - 1 ;z) - (1 - 
z 2 )F(z~';z). From this, we also find expressions using incomplete 
elliptic integrals for equations <A3t and <A4> : 



S[z c ] 



Irl 2 



' r E (> 



S[ «=2f(1 ;W )-4e(1 ;M 
for equations <A8t and <A9t ; 



S[z m ] 

dS [z,„ 
dm 



Irl 2 



El'—- — 
~ 1 1 |y|' m 

2m 2 



Irl 2 



Irl 2 



El ™;JZl) +F ( — ;^ 

m 2 \ |y| ' m J \ |y| ' m j 



and for equation JA1 1> 



S[z p 



Irl 2 - 1 



FI M ;M 



E 5f w 



(A12) 
(A 13) 

(A14) 
(A 15) 

(A16) 



IPress et ail ll992t) incorporated the so- called ICarlsorl 1 19771) 
symmetric form of the elliptic integrals, which may be advanta- 
geous for certain considerations over the conventional Legendre- 
Jacobi standard form. Two of Carlson's integrals are used to re- 
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express the elliptic integrals of the first and the second kinds: 
R F (x,y,z) = ijT (t + x)- >l2 (t + y)- i/2 (t + zr i,2 dt 

R D (.x,y,z) = | J (t + x)- l/2 (t + y)- l/2 (t + zr 3l2 dt. 

From the Jacobi form of the integrals given above, by changing 
the integration variable to / = (s/v) 2 - 1, it is easy to establish the 
transformations: 

F(s;k) = sR F (l-s 2 ,\-k 2 s 2 ,\) 

EO;ifc) = sR F (l - s 2 , 1 -k 2 s 2 , 1)- — R D (l - s 2 , 1 -k 2 s 2 , 1), 

and subsequently, K(k) = R F (0, 1 - k 2 , 1) and E(k) = R F (fl, 1 - 
k 2 , 1) - (k 2 /3)R D (0, 1 - k 2 , 1). Using this, we can re-express equa- 
tions (KT\ and (A2\ 

5fc] = T q^[^(0,c 2 ,l)-^ D (Q,c 2 ,l) 

5 [fij = 2ji + ^lrl 2 ^ (0, c 2 , l) - 4tf„ (0, c 2 , l) 
where c 2 = 1 - |y| 2 , equations JA3> and JA4> 

^ fe] = ,,.,„..!, „ (0. g 2 , i) - (o, g 2 , l) 



IrKIrl 2 - 1) 



l 

M 



-R D (0,q 2 ,l)-2R F (0,<r,l) 



where q 2 = 1 - \y\ 2 , equations <A5> and JA6> 

2m 



m 2 - |y| 2 



dm (m 2 - |y] 2 ) 2 



^(0.4.1) -^(o.4.i) 



drl 2 * D (o, 4. i) - (™ 2 + Irl 2 )^ (o, 4- 1) 



where c 2 = 1 - (|y|/m) 2 , equations <A8> and <A9t 



* F (0,4,l)--*i,(0,«£,l) 



SUJ IrKIrl 2 -™ 2 ) 

dS [z m ] m 



dm lyldyl 2 - m 2 ) 2 
where q 2 m = 1 - (m/|y|) 2 , and equation <A1 1 

S[z p ] 1 



(|y| 2 + m 2 )R F (0, <7 2 , l) - -m 2 R D (0, 9 2 , l) 



3|yl(lyl 2 - 1) 



tf D (o,<7 2 ,i). 



APPENDIX B: THE NUMBER OF IMAGES OF A 
RATIONAL LENS MAPPING 

In this Appendix, we examine the problem: what is the maximum 
number of images for N point masses with shear? In fact, it is 
scarcely any more work to consider the more general case of the 
deflection function that is a rational funct ion of complex posi- 
tion. T his problem was in fact solved by iKhavinson & Neumannl 
(2006) under still more general considerations. In particular, for 
the case of N point masses in the absence of shear, their result 
provides us with an affirmative answer to the conjecture by Rhie 
(2001, see also Rhie 2003): the maximum number of images for 
N point masses ( N > 2) with no shear is 5(N - 1). However, 
while the result of Khavinson & Neumannl i2006t) is strictly true, 
the limit can be further lowered dependin g on the particular be- 
haviou r of the rational function, as seen by Khavinson & Swiatek 
1 2003). Here, we revisit the proofs of Khavinson & Swiatek 1 2003) 



and Khavin son & Neumannl 120061) . and derive specific limits de- 
pending on the asymptotic behaviour of the rational lens mapping. 
Let us suppose that the lens equation is written as 

£ = z - s(z) 

where s(z) is a rational function of degree r > 2 - i.e, a rational 
function that is neither a constant function nor a Mobius transfor- 
mation. 5 Therefore, s(z) '■ — > C m is surjective - 'onto' - but 
not injective - 'one-to-one'. Here Coo = C U jooj is the 'Riemann 
sphere.' By considering the complex conjugate, the lens mapping 
can also be written, after rearranging terms, as 

z = f(z) = 1+ s(z) 

where f(z) is also a rational function of degree r, and in addition 
f'(z) = s'(z)- Then, the image positions are found among the fixed 
points of a degree-r 2 rational function h(z) = (f°f)(z) = f[f(z)] or 
equivalently the zeros of a rational function g(z) = z - h(z). Here, 
f(w) = f(w) is also a degree-r rational function of its argument so 
that h(z) is in fact a degree-r 2 rational function of z. On the other 
hand, the degree of g(z) is either r 2 + 1 if lim-^„ h(z) is finite or 
r 2 if it is infinite. Actually, if lim^,*, s{z) = °o, then h(z) diverges 
as z — > oo, whereas linu^c h(z) is finite if lim^^, s(z) = a is finite 
except for the case when s(z) has a pole at z = <f + a. Here, the 
incident of the exceptional case explicitly depends on the specific 
choice of the source position (£), and so this cannot be a generic 
occurrence for any lens system with an arbitrary source location. 
Hence, ignoring the exceptional cases, we find that the degree of 
g(z) is either r 2 if s(z) diverges as z — > °° or r 2 + 1 if otherwise. 

Now, according to the theorem due to Pierre J. L. Fatou (1878- 
1929) and Gaston M. Julia (1893-1978), we have 

Theorem B.l. (Fatou-Julia) "Every attracting cycle for a ratio- 
nal function (degree r > 2) attracts at least one critical point," or 
equivalently "Every basin of attraction for a rational map (degree 
r > 2) contains at least one critical point." 

The proof requires some familiarity with complex analysis, in par- 
ticular, Schwarz's Lemma and/or the Riemann Mapping Theorem, 
which is beyond our scope (see e.g. lCarleson & Gamelinl 1993h . In- 
stead, here we simply accept the implications of Theorem lB.ll for 
the current problem as a fact without detailing any proof: For any 
attracting 6 fixed point of h(z), there is at least one critical point 7 
('stationary point' or 'ramification point') of h(z) that converges to 
the fixed point under the iterative mapping of h(z), that is to say, 

Theorem B.2. For any point zt 6 C„ that satisfies h(zf) = Zt and 
\h'(zt)\ < 1, there is at least one point z c 6 C„ such that h'(z c ) = 
and lim„^„o h"(Zc) = Zt. 

Here, h 2 (z) = {ho h)(z), h 3 (z) = (h o h 2 )(z), and so on. We need 
to take particular care in the interpretation of h'(z) when z = °° 
or h(z) = oo. By considering its composition with the map z~', we 
deduce the following: 



Formally, this indicates that the image .y(C TC ) of the Riemann sphere under the given (meromor- 
phic) function s(z) : Ceo — > Ceo is anr-fold covering of Co,. That is, the number of solutions in C TO of 
s(z) = a (counting multiplicities) is r for any fieC. With a rational function s{z) = ^n(-)/-^d(^) where 
s n (z) and sjiz) are polynomials with no common factor, this is basically /' = max{deg ,9 n (z),deg s^iz)} 
where deg p(z) is the degree or the order of the polynomial piz). 

6 The fixed point of h(z) is said to be super-attracting, attracting, neutral, or repelling, if h'(z) = 0, 
\h'(z)\ < 1, = 1, or |/i'(z)| > 1, respectively. 

7 

The critical point of a function is the point at which the mapping given by the function is locally 
stationary. To avoid any confusion with the critical points of the lens mapping, we will avoid the use 
of this term in favor of the stationary point. For an analytic function s{z), they are equivalent to zeros 
of s'(z). However, for a meromorphic function s(z) of Riemann sphere onto itself, this identification 
should be applied with a caveat due to the association of = 0. That is, the point at which 
s'(z) = with n is a positive integer is not necessarily a critical point. For example, the mapping 
s(z) = z -1 is regular at z = 00 although s'(z) - -z~ 2 . In fact, for the discussion in this appendix, the 
critical point may be better understood as being a ramification point in the topological sense. 



© 2006 RAS, MNRAS 00Q.fTWT9l 



16 An & Evans 



(i) If h(oo) = lim z ^ 00 A(z) = oo, then /i'(oo) = 
lim^oo z 2 h'(z)/[h(z)] 2 . Specifically, if h(z) ~ az as z — » oo with 
a + being a constant, then /i(oo) = oo and |/i'(co)| = |a| -1 , and if 
h(z) ~ az" (n > 2) as z — » °°, then r?(oo) = oo and /i'(oo) = 0. 

(ii) If |/i(oo)| = | lim.^oo h(z)\ < oo is finite, then h'(oo) = if and 
only if lim.^„ z 2 h'(z) = 0. 

(iii) If h(z p ) = oo for a finite z p , then h'(z p ) = if and only if 
lim z _, Zp h'(z)/[h(z)] 2 = 0. In other words, if z p is a higher-order pole 
of h{z) than a simple pole, then h'{z p ) = 0. 

In the following, we want to prove that Theorem lB.2l implies that 

Theorem B.3. For any point zo e C m that satisfies f(zo) = Zo and 
\ s '(zo)\ < 1» there is at least one point z c 6 Coo such that s'(Zc) = 
and lim„^oo(/ ° f)"(z c ) = Zo- 

To prove this, we first observe that h(zo) = (f ° f)(zo) = flf(zo)] = 
f(zo) = /(Zo) = Zo, that is, zo is a fixed point of h(z). Next, from the 
chain rule for the derivative of a composite function, 

h\z) = (fof)'(z) = f'[f(z)]f'(z) = s'[f(z)]s'(z), (Bl) 

we find that h' (zo) = S'[/(zo)Mzo) = s'(zo)s'(zo) = s'(z )s'(z ) = 
l- s ''(zo)l 2 and so that |/i'(zo)l < 1- (Here, s(z) = s(z) or equivalently 
s(z) = f(z) - <T-) Then, from Theorem IB. 21 for each zo, there ex- 
ists at least one point z c such that h'(z c ) = and lim„^ m h"(z c ) = 
lim,,^oo(/ ° f)"(z c ) = Zo- From equation IB 1 1 , the h'(Zc) = indi- 
cates that we have either s'(w c ) = where w c = /(z c ), or s'(z c ) = 0. 
If the latter is the case, z c is simply the point that we are looking 
for. If the former, we find that s'(w c ) = s'(w c ) = and that 

lim(/o/)"(w c ) = lim(/o/)»[7fc]] = limC/o /)"[/&)] 

n— >oo n—>oo n^oo 

= lim /[(/ o /)"(!,)] = / [lim(/ o /)"(z c )| 
= /[lim(/o/)«( Zc )] = f(zo) = J(^) = zo. 

Here, we have exploited the associativity of the composition of 
functions and the commutativity between the limit and the any ana- 
lytic map. We conclude that w c is the point that we are looking for. 
While it is not of direct interest here, following similar arguments, 
we can further show that, for any point zo satisfying the condition 
part of Theorem lB.3l there are at least (r + 1) points z c - where r is 
the degree of s(z) - such that h'(z c ) = and linv-«o h"(zc) = Zo (i-e., 
the consequence part of Theorem lB.2> . In particular, at least one of 
them satisfies s'(z c ) = 0, while r of them map to a single point vv c 
under / - i.e., they are all of the r pre-images of w c - where vv c 
satisfies s'(w c ) = 0. 

Note that the condition part of Theorem lB.3l is actually equiv- 
alent to zo being an actual image position with ^f(zo) > 0. That is, 
f(zo) = I + s(zo) = Zo so that f = zo - Kzo), which is just the lens 
equation, whereas J (zo) = l-k'(z )| 2 = (1 - k'(zo)l)(l + W(zo)\) > 
if (and only if) |s'(zo)l < 1- Hence, the immediate corollary of 
Theorem lB.3l is that 

Corollary B.4. (Rhie 2001) The number of positive parity images 
(n + ) for the lens mapping £ = z — s(z) that can be described by a 
rational function s(z) of degree r ^ lis bounded by the number of 
stationary points of s(z). 

However, it is known that the latter number is exactly 2r - 2, if the 
point at infinity is included and the multiplicity has been counted, 8 

8 This follows from the Riemann-Hurwitz formula. That is, since s(z) : Ceo — » C TO is an analytic 
map of degree r of the Riemann sphere (whose Euler characteristic is 2) onto itself, the total number 
of ramification points (or stationary points) is 2r - 2. In particular, for a polynomial p(z) of order r, 
its derivative is a polynomial of order (r - 1) so there are (r - 1) stationary points in C. In addition, 
since p(z) ~ az r (fl ^ 0), the point at infinity is degenerate (r - 1) times, and thus it contribute an 
additional (r - 1 ) counts of ramification points. 



and therefore, n + < 2{r - 1) <Rhidl200lHKhavinson & Neumann! 
2006). 

While this limit is strictly correct, the number of distinct sta- 
tionary points of s{z) may be smaller than 2r - 2, if there are de- 
generate solutions for s'(z) = 0. In particular, if there exists a point 
Zc such that ,s'(z c ) = . . . = s ( "- [) (z c ) = and .v <n) (Zc) * or if s(z) 
possesses an n-th order pole, with n > 3, the number of distinct sta- 
tionary points and consequently the limit on is lowered by n - 2 
for each such point. Another possibility for lowering this limit fur- 
ther depends on the asymptotic behaviour of s(z), provided that we 
are only interested in images with |zol < oo. 

• First, if lirn^oo s(z) = a is finite, z = oo is not a fixed point of 
h(z) provided that s(z) does not have a pole at £ + a. On the other 
hand, if s(z) ~ a + bz~" where b t is a constant and n > 3 is an 
integer, then s'(oo) = ... = .^"-"(oo) = and ^"'(oo) t so that 
z = oo is an (n- l)-times degenerate stationary point. Consequently, 
we lose (n - 2) distinct points that might have iteratively converged 
to one of the positive parity images, and thus, the upper limit on n + 
is lowered to n + < 2r - n. 

• Secondly, if s{z) diverges at least quadratically, we find not 
only that .y'(oo) = h'(oo) = 0, but also that h(oo) = oo. In other 
words, z = oo is a stationary point and a (super-attracting) fixed 
point. For this case, the point z = oo obviously does not itera- 
tively converge to any finite fixed point, and therefore, this point 
should not be considered for the limit for n + at all (counting all 
of its multiplicity). If s{z) ~ az" (a t and n > 2), we find that 
i-'(oo) = . . . = / n-1 >(oo) = and 5 (n) + 0, and thus, the limit is 
lowered to n+ < 2r - n — 1 . 

• Finally, if s(z) diverges linearly so that s(z) ~ az (a t 0), 
then .y'(oo) = aT l + and a'(oo) = /(oo) = oo and consequently, 
h'(oo) = |a|~ 2 and h(oo) = oo. With the identification of oo = oo, 
we can conclude that there is at least one finite point such that 
j'(Zc) = and lim„^„ h"(z c ) = oo if \a\ > 1 since z = oo is a 
non-stationary attracting fixed point. Therefore, for \a\ > 1, there 
is one less point with s'(z) = that can iteratively converge to one 
of the finite positive parity images, and consequently, the limit is 
given by « + < 2r - 3. However, if \a\ < 1, the fixed point at infinity 
is repelling and thus does not affect the behaviour of the iteration 
of finite points, and the original limit n+ < 2{r - 1) is maintained. 

Bl Index Theorem for a rational lens mapping 

The limit on n + can be extended to the limit on the number of total 
images (n + + nJ) by deriving the relation between the numbers of 
positive and negative parity images (n+ and respectively). Let 
us first think of a function F(z, z) = s(z) — Z — X and 

J( „ 9-Fdz+9=Fdz i'(z)dz-dz 

d In F = = : — - . 

F s(z) - z - f 

Then, the contour integral 

JT s'(z)dz - dz 
iA flc argF=(p — — z — - 
Jdc s(z) - z - < 

along the boundary of a domain C c C vanishes, provided that the 
domain C does not contain any image position (i.e., zeros of F) 
nor poles of s(z) [i.e., s(z) t z + I and \s'(z)\ < oo for v z 6 C], 
This follows because the integrand is exact, so that it is closed in C. 
Furthermore, for an arbitrary domain C that may contain image po- 
sitions and/or poles, the integral can be evaluated by collapsing the 
contour to infinitesimally small circles around each image position 
and pole. 

Suppose that zo is an image position that is not on the critical 
curve. Then, s(z) - I + Zo + (z - ZoK(zo) + &(\z - Zol 2 ) and s'(z) - 
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Table Bl. The number of non-infinite regular images for a rational lens mapping of degree r > 2. Here, the lens equation is given by ( = z. - Hz) with s(z) 
being a rational function of z, of degree r > 2. The asymptotic behaviour of s(z) as z — > °° determines the maximum possible number of non-infinite positive 
parity images. The positive integer n is basically the ramification index at infinity. If s(z) is finite at infinity, the total count of the order of poles in C is the same 
as the degree r while it is r - n if s(z) diverges at infinity. Physically, the number of poles corresponds to the number of point masses. The positive integer d is 
the duplication index, which counts the number of ml> 3)-th order poles as well as points such that s'(z c ) = ■ ■ ■ = s 'fc) = and s' m '(Zc) + with m > 3. 
For each such point contributes m - 2 to d. 



Asymptotic Behaviour 
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s'iza) + G{\z — Zol), and thus, the integral along the small circle 
z = Zo + ee 1,p (0 < ip < 2n) around zo reduces to 



iA, 3C argF = i 



r 2n d s'(z y* + e-'" + ■ 
Jo V s'(zo)^ - e-* + ' 



(£) 
(0' 



and can be evaluated as the limiting value for e — » 0. That is, if 

k'(zo)l < 1, 



f 2 " 5'( Zo )e^ + e- 
Jo 



s'(z )e'f - e-'f 
whereas if |.s'(zo)l > 1, 

-271 



Jo ^ 1 - ,,(z )e 2 - = ~ 2n > 



C- n d s'(zo)^ + e-^ p d 1 + [Azo)T 
J V ^'(zo)e^ - e-V J V 1 - [5'( Z0 )]- 



: 27T. 



Here, both results can be obtained by series-expanding the denom- 
inator and using J Q e'" v di/) = where n is a non-zero integer. 

Similarly, for a pole z p of s(z) of order n, by factoring out 
the diverging part, it is possible to write s(z) = (z - z p y n s v (z) 
where s p (z p ) is non-zero finite and s'(z) = -n(z - z p r"~' s p (z) + (z - 
Zp)~"s p (z). Consequently, the integral along the small circle around 
the pole z = z p + ee v (0 < y < 20) is given as 



i A ac arg F = i 



271 _ -ns p (z) + ec^s' p (z) + e" +l e ( "- ly "f 
i p (z) - e"e"^(| + z p ) - 6" +1 e<"-"^ 



d<,p ■ 



the limiting value of which as e — > is -2«7ti, assuming 
lim z _, Zp i p (z) is finite. Hence, for any domain C chosen large enough 
to contain all of image positions and poles of s(z) in C, we find that 



Age arg F = 27t(n_ — «+ - n») 



(B2) 



where n + and w_ are the total numbers of positive [|j'(zo)I < 1 so 
that ^(zo) > 1] and negative [|.?'(zo)l > 1 so that ^(zo) < 1] parity 
images, respectively, and is the sum of orders of all poles. This 
result is in fact a particular case of the argument principle applied 
on the function F. 

However, for a large enough contour dC, the integral can also 
be evaluated using a contour around the point at infinity 



iA flc argF 



.s'(z)dz - dz 
s(z) -I- 1 



w s'(w )dw — w dw 
s(w- l )-w->-Z 



The limiting result as w — > is dependent on the asymptotic be- 
haviour of s(z) as z — » oo. If lim ; ^„ s(z) is finite, using the asymp- 
totic expansion form s(z) ~ a+bz~" where a and b t are constants 
(a may be nil) and n is a positive integer, the integral reduces to 



r dw r- n 

i Aac arg F = (J) — = -t dtp = -2m. 

Joe w Jo 



JdC w JO 

On the other hand, if s(z) diverges at least quadratically, that is 



s(z) ~ az" with n ^ 2 and a # 0, then 
dw 

he w Jo 
Finally, if s(z) ~ az(a t 0) diverges linearly, then 



i A flc arg F = n 



r dw _ r 2n 

Jac w Jo 



dtp = 2n7ti. 



i A, 3C arg F 



as -up + e up 
dip : — 

o ae,- l f - e 1 * 



2m if \a\ > 1 
-2ni if \a\ < 1 



By summarizing the results, we arrive at a relation for the difference 
between the numbers of positive and negative parity images 



n m - I if lim z ^„ \z 1 s(z)\ < 1 
n„ + n = r if lim^oo |z _1 i(z)| > 1 



(B3) 



where n is the ramification index of ,v(z) at z = oo. Since s(z) is 
a rational function, if s(z) diverges as z — > oo, n is the same as 
the asymptotic power index of s(z), which must be an integer for 
this case. Furthermore, n m is in fact the same as the order of the 
polynomial in the denominator when s(z) is expressed as a quotient 
of relative prime polynomials. Since the asymptotic power index 
for a rational function is basically the difference of the orders of 
those polynomials, and the degree of the rational function is the 
larger of two, we find that r = n m + n if ,v(z) is divergent as z — > oo. 
This result is in fact the extension of the so-called index theorem 
( Burke fl98lTl to the case when the deflection function is given by a 
rational function. 

The fact that n_ - n+ is a fixed number for given s(z), in par- 
ticular, being independent of f, further implies that for a given lens 
system, change of source position can only create or destroy images 
of opposite parity in a pair. The upper limit on the total number of 
non-infinite images can be directly found from the limit on n + using 
n + +n_ = 2n + + (n_ - n + ). The final result on the limits of number of 
images for a number of rational lens maps are summarized in Table 

IbTI 

Comparing the upper limit on the number of images to the 
degree of the imaging equation g(z) = z - h(z) = 0, which is either 
r 2 if s(z) diverges as z — > oo or r 2 + 1 if otherwise, we find that they 
can be same only if r = 2 or r = 3 [i.e., r 2 - 5r + 6 = (r - 2){r — 
3) = 0]. If s{z) is a non-degenerate Mobius transformation, that 
is, a degree- 1 rational function, explicit calculations can establish 
that the imaging equation reduces to either a linear equation if j(z) 
is a linear polynomial, or a quadratic equation if otherwise, with 
all of its solutions yielding image locations. In other words, the 
imaging equation must have solutions that are not image locations 
if r > 4. This also indicates that the moment sum invariants that 
may be derived through the similar methods as in section l3~T1 will 
be valid only if the degree of the rational deflection function is 1, 
2, or 3, since, if otherwise, there must be spurious roots. 
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APPENDIX C: THE MEAN MAGNIFICATION FOR THE SOURCE INSIDE THE CAUSTICS 



As noted in section l5~2l the ratio of the areas under the caustics and the pre-caustics gives the mean magnification for the source inside the 
caustics. So far, we have derived the areas under the caustics (eqs, !55l and lA2> for < |y| < 1 and (eqs. l60lfA4l and lA13> for |y| > 1, and under 
the pre-caustic (eqs, l86llAlll and lA16> for |y| > 1. To complete the discussion for all values of |y|, we need to figure out the area between the 
two pre-caustics for < y < 1 (eqs. l63l and !64> . As it turns out, this can be expressed in terms of the standard elliptic integrals (all four forms 
are equivalent): 

4|y| 

S [z p + ] - S [z p ] = [(1 + cJILQi, x) + K(x) - 2cE(x)] 



4|rl 

c 5/2 



lrl 2 d + 2c) 



l+c 

4|y| 
(1 + 2c) 1 
4|y| 



Tl(-A, x) + c(l + 2c)K(x) - 2cE(x) 



c2(1 ™ lj2 [4|y| 2 n(y,y) + 2c(l + 2c)K(y) - (1 + 2c)E(y)] 
m [4cH(-t, y) + 2K(y) - (1 + 2c)E(y)] , 



c 2 (l +2 C y 

or in terms of the Carlson symmetric form: 
S lz + J - S [zT] = % ((2 - c)i? F fo, -, l) - \ 



4|y| 

r 5/2 



^Jo.-.l.-U 



4|y| 2 



(1 + 2c) 



tf c 0,-,l 



/? F |0,-,1| + - 



(l+2c)|y| 2 



(1 + c) 2 



0,-,l,tf + 



3 - 4|y| 2 
(l+2c)' 



r d o,-,i 



c 2 (l + 2c) 1 ' 2 

Here c = ^1 - Irl 2 , 
1 



4|y| 



2(1 + c) 



Irl 2 



p = — = 1 - T 



2c(l + c) 
1 2c 



3 - 4|y| 2 
4c(l + 2c) 



Irl 2 



l 



4|y| 2 



8= - = 1 -A 



1 



(1 + c) 2 
1 +2c 



1 +2c 



(l+2c) 2 



1+yU 1+c' v 1 + v 2(1 + c) ' 

and H(n, A:) is the complete elliptic integral of the third kind: 

.71/2 

n(n,k) = 



, 6 1 + 2c 

1 - x = - = : 

p Ac 



. 2 p 4c 

1 - y = - = 1 

<5 1 + 2c 



(l +nsin 2 (f>) yjl -k 2 sin 2 (p 
and Rj(x, y, z, p) is an additional Carlson's integral: 

Rj(x,y,z,p) s I f (t + x)- l l 2 (t + y)- ll2 (t + zy U2 (t + p)- l dt. 
2 Jo 

Note that Tl(n,k) = R F (0, 1 - k 2 , 1) - (n/3)Rj(0, 1 - k 2 , 1, 1 + n). With these results, we can find the analytic expressions for the mean 
magnification for the source inside the caustics. 

As a simple illustration of the use of these expressions, we derive the asymptotic expansions of (M)/M^ in , which are also plotted in the 
lower panel of Fig. [4] 



(M) 
M c . 

mm 

w 

min 



(M) 

min 



(M) 
M c . 

min 



3tt 
1 - 



3K - -2E - 



bf 
In 



c 

+ — 



13K(I)-14E(- 

37T-8 



Irl 4 



288tt 
10- 3tt 



123K - -98E - 



4-Tt + 21n(c/4) 2 ( 4 - n + 2 ln(c/4) [4 - n + 21n(c/4)] 2 



+ —<- + 



157T-41 



267tt - 469 - 36tt 2 



(10-37T) 2 



8 [2 2[4-7t + 21n(c/4)] 4[4 - tt + 2 ln(c/4)] 2 [4 - n + 21n(c/4)] 3 



+ ^(lyl 6 ), (Irl -» 0) 



(c 6 ), (|y| -> I") 



1 - 



1 



2 + \n(q/4) 4 



--10- 



14 



2 + lnfo/4) [2 + ln( ? /4)] 2 



392- 



724 



357 



72 



2 + 



128 
1 

4~jyp ~ 64|y| 4 



1 



2 + \n(q/4) [2 + ln(g/4)] 2 [2 + \n(q/4)V 



5 ) (|y| -> oo). 



+ 0(q\ (|y| -> 1 + ) 
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Here, c 2 = 1 - \y\ 2 and q 2 = 1 - |y|~ 2 . For the first two cases (0 < \y\ < 1), we have M c . = [|y|(l - |y| 2 )] -1 while for the last two 

drl > i), = (2\ 7 \ 2 - i)/(| 7 | 2 - i). 

This paper has been typeset from a TpX/ ETgX file prepared by the author. 
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